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Abstract 



We examine the static potential in the presence of flavors in the perturbative back- 
reacted D4/D8 system from localized D8 branes, focusing in particular on the Sakai- 
Sugimoto model. For the case of long strings we find the flavor corrections to the static 
potential which are of exponential form. We then investigate shorter Wilson loops and 
express their energy analytically in terms of the lengths of two neighboring Wilson loops. 
Moreover, we find that beyond a certain scale the static force in the backreacted back- 
ground is reduced compared to one in the probe limit, as expected due to screening effects. 
We also compare the string world-sheets in the two backgrounds and find how they get 
modified by the backreaction. Our results are supported by numerical computations as 
well. Finally we discuss our results in comparison with the lattice data and comment on 
the issue of physical scales which seem to lie at the heart of the reason that obstructs our 
model at this level to fully describe QCD. 
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1 Introduction 



The Wilson Loop (WL) is a an object that contains fundamental information about a 
gauge theory at a non-perturbative level. Its ground state gives us the static potential, 
i.e. the interaction between two infinitely heavy quarks. For SU(N C ), the pure gauge 
theory at zero temperature has a unique, confined phase with a static potential of the 
form 

c gauge 

E = const. + ah + — 1 (1) 

Lj 

Each of the above terms needs a separate discussion. The constant is non-physical. Typ- 
ically one considers the static force dE/dL to get rid of it. The string tension a for 
small N c , can be computed in lattice Monte Carlo simulations as the slope of the static 
potential, as a function of the static quark separation L. It has dimension of length -2 
and for that reason it is usually expressed in terms of the lattice spacing or better in 
terms of some other physical dimensionful scale. Of special importance is the Coulomb 
like term because its coefficient is dimensionless. In the ultraviolet (UV), that is at small 
L, it is related to the asymptotically free coupling of the gauge theory and evolves ac- 
cording to the well known, perturbative renormalization group equations. Towards the 
infrared, it stops running roughly around a particular scale, the "Sommer scale" and 
from there on it assumes a constant value. This constant value has been predicted to be 
cf auge = — (d — 2)7r/24 in ^-dimensions and for any N c [1]. The prediction derives from a 
picture which has a non-critical bosonic string connecting the two static quarks and the 
(Gaussian) fluctuations of the associated world-sheet in a flat transverse space, give rise 
to this "Luscher term". Even more impressively the value was successfully confirmed in 
a Monte Carlo simulation [2]. Finally the dots in Eq. (1) stand for unknown corrections, 
at least from the point of view of the lattice. 

Things get more involved when dynamical quarks are added to the gauge theory. To 
begin, the string tension can be measured straightforwardly via Monte Carlo simulations, 
as in the pure gauge case. On the other hand, to our knowledge there is no string derived 
prediction for the flavor correction to the Luscher term up to date, at least for QCD. 
As far as the lattice is concerned the only estimate of a correction to the QCD static 
potential which could perhaps support the effort to extract such a prediction from string 
theory is a WL analysis with Nf = 2 dynamical flavors that appeared recently [3]. The 
result of this study is that the value of the Coulomb term coefficient C\ is consistent with 
the prediction of the Cornell model [4] c\ ~ —0.52, which means that the Nf = 2 flavors 
approximately double the pure gauge value. What is not clear at all is if this correction, 
which we call from now on cf , is purely geometrical like it is in the pure gauge theory or 
dynamical. From the field theory point of view nothing ensures that a stringy description 
even exists in QCD in the presence of dynamical flavors. We obviously suspect that there 
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may be such a description and we study it in this paper. 

Regarding the higher order corrections in Eq. (1), they are even harder to probe in 
lattice QCD unless perhaps an independent input, for what kind of terms one should 
expect is found. String theory offers us the possibility to compute such corrections, at 
least to the pure gauge static potential. Already in flat space, an infinite number of terms 
contribute corrections that come as inverse powers of L. The form of this polynomial series 
depends on the restrictions that one puts on the derivative expansion of the world-sheet 
action. The most general expansion is obtained when only Lorentz invariance constrains 
the world-sheet action. If on the top of it one assumes modular invariance, one obtains a 
restriction of the above series. 

The question that naturally arises is which is the string that corresponds to the QCD 
WL. The exact answer to this question is not known but we do know a class of espe- 
cially promising approximations via special, curved holographic backgrounds using the 
AdS/CFT correspondence [5,6]. The first such computation was that of the static po- 
tential in the Af = 4 SU (N c ) SYM, in the limit where A = gy M N c is large and fixed and 
N c — y oo. In this setup, the WL is the boundary of a world-sheet that extends in the 
bulk of a curved (AdS 5 ) space-time. The static potential was derived from the minimiza- 
tion of the world-sheet surface with the appropriate boundary conditions and yielded a 
Coulomb-like potential [7, 8] , as expected from a conformally invariant gauge theory. 

By now there is a good understanding of many extensions of the original Maldacena's 
proposal. The AdS/CFT correspondence has been extended to theories with less or zero 
supersymmetries, as the (3 deformed ones 1 and the quiver theories with the Sasaki-Einstein 
gravity duals. On the other hand there have been found also several confining backgrounds 
at zero and finite temperature, aiming to provide the ground to understand better QCD. 
One important development in this direction was the understanding that fields which 
transform in the fundamental representation as the QCD quarks do, can be introduced 
by adding appropriate flavor branes in the gravity dual background [9] . Usually this is 
done in the probe limit where the backreaction of the flavor branes is neglected, and by 
studying the worldvolume of these branes one can observe in the gravity dual background 
many physically expected phenomena and to compute for example the meson spectrum 2 . 
Nevertheless, in the quenched approximation the mesons are affected from the glueballs 
but not the other way around. So ideally one would like to consider the backreaction 
of the flavor branes to the original background, in order to obtain more realistic models. 
The equations needed to be solved for the derivation of the backreacted backgrounds 
are in general second order, non-linear partial differential equations, which include Dirac 

x The Wilson loop for the deformed theories at least for large N, is expected to be similar to one of 
the TV = 4 SYM theory [10]. 

2 See [11] for a review on the AdS/CFT with flavors. 
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delta functions due to the localization of the flavor branes. To solve these equations it 
is a challenging problem. On the other hand one can consider an approximate approach, 
and allow a smearing of the flavor branes instead of a localization. This approximation 
simplifies the differential equations, but the smearing breaks a U(Nf) flavor group to 
U(l) Nf and hence theories with U{Nf) flavor groups can not be examined exactly 3 . 

In this paper, motivated from the previous discussion about the lattice approach and 
the difficulties encountered in AdS/CFT away from the probe limit, we study how the 
static heavy QQ potential is modified in holographic QCD when the backreaction of the 
flavor branes is taken into account. For this reason we use the Sakai-Sugimoto [13] model 
which is the closest one so far to the holographic description of large N c QCD. Taking 
advantage of the fact that the backreaction of the localized branes has been found near 
the tip of the geometry for the D8 and D8 branes placed on the antipodal points of the 
compactified circle [14] , we argue that the leading corrections to the static quark potential 
for the particular WL we consider come from this region. Doing that we achieve to study 
analytically as well as numerically the flavor corrections on the static potential and the 
static force. The corrections to the potential in the non-flavored background for very 
long strings are expected to be exponential [15], and we show that the flavor corrections 
in this limit are also exponential although with different factors. As a result, for such 
strings, we do not find any correction to the string tension or to the Luscher term. In 
order to compute the static potential of short strings we perform a series expansion which 
converges to a string length of order one. We compute the static potential there and 
we find flavor corrections to the string tension but no corrections in inverse powers in 
the string length, beyond those that come from the a' expansion. In several cases we 
compute explicitly all the interesting coefficients appearing in the static potential. We 
perform a first comparison of the static force in the holographic approach with lattice QCD 
data. A qualitative similarity in their behavior is observed but as we will demonstrate, a 
quantitative comparative study is subtle. 

The paper is organized as follows. In section 2 we review the Witten-Sakai-Sugimoto 
background (WSS) and the localized backreacted solution we are interested in. In section 
3, we first derive the flavor corrections to the energy in the limit of large distance between 
the QQ and find them to be exponential as in the probe approximation but with different 
factors. Then we present a method which allows us to compute short an intermediate 
distance corrections. The numerical analysis is included in section 4 where we explain 
how we perform it in order to be valid for our approximation and we present several 
results for the static potential, the color force and the string worldsheet comparison in 
the two backgrounds. Finally, we summarize and discuss our results in the last section 

3 See [12] for a review on the backreaction of smeared flavored branes. 
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and include two appendices to support the analysis of the main text. 



2 The general setup 

2.1 The Witten background for D = 4 

The unflavored version of the WSS model, is similar to the background derived by Witten 
[16] and consists of N c D4 color branes, where one spatial direction is compactified on 
S* 1 with radius p. The submanifold spanned by the radial coordinates u and has a 
topology of a cigar which ends at u — Uk, where to avoid the singularity at the tip, one 
requires that 

x 4 ~x 4 + 2irp, p : = (2) 
The background dual to the 5-d pure Yang- Mills theory is the D4 brane metric solution [16] 
ds 2 = (|) ^ (i^dzW + f(u)dxl) + (|) 3/2 (J^ + u 2 dn^j , (3) 

= g s Q) 3/ \ ^ = dC, = ^ e 4 , f(u) = 1 - ^ , R 3 = 7cg s NJ 3 s , 

where /i — 0, .., 3 and (ifi 2 , 64 and V4 = 87r 2 /3 are the line element, the volume form and 
the volume of a unit S 4 , respectively. Also, the string length is l s = y/a'. The couplings, 
satisfy the following relations: 

Clearly, when p — > 00 we reach the (UV) limit where the dual gauge theory is 4+1 
dimensional and in the opposite (IR) limit we expect to see an effective 3+1 dimensional 
theory, possibly with some memory of its five-dimensional origin. It can be shown [17, 18] 
that in the IR the static potential, consistently with dimensional reduction, develops a 
string tension equal to 

- (0) = ^7f^) 3/2 - (4) 



2na' \ RJ 

We will repeat and extend the proof of this later. It was found that when a fair number 
of higher derivative terms are taken into account in the world-sheet action, the Luscher 
term retains its universality but the string tension does get modified [19]. We review the 
results regarding a' corrections in Appendix B. 

In the IR limit the static potential is found to be 

E m = const. + a^L - ^ + 0{L^) , (5) 
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where is the renormalized string tension and the Coulomb like correction is the leading 
correction in an a' expansion, the above mentioned Liischer term. 

The computation of the static potential in the UV is left for section 3, along with 
a possible alternative approach towards a comparison with lattice gauge theory. From 
Eq. (40) we find a 5-d Coulomb potential. Thus we arrive at a picture of the phase dia- 
gram of the Witten model where the 5-dimensional SU (N c ) gauge theory passes, without 
crossing a sharp phase transition, from a Coulomb phase to a 4-dimensional confined 
phase, as the theory gets dimensionally reduced when p — > 0. The theory never looses its 
5-dimensional memory since all light states are proportional to 1/p and p can never really 
be zero. In addition, as noticed in [20], it does not loose its supersymmetric origin either, 
since there are fermionic mesons originating from the world-sheet fermions that are also 
set by the same scale. We will not compute here flavor corrections in the UV limit, since 
their immediate importance is questionable. 

The puzzle that arises then from all this is that the correction to the Coulomb term 
that the lattice sees [3] and which is of order c\ ~ 0{— 1/4) can not be explained by 
Witten's background even if we include a' corrections. Moreover, adding probe flavor D8 
branes can not improve the situation because in this limit the background, and therefore 
the WL, is unaffected. The only left over possibility we see is to consider the backreaction 
of the flavor branes and recompute the WL in the corrected background. 

2.2 The WSS background with backreaction 

The Witten-Sakai-Sugimoto model is the extension of the model of the previous section 
with a U(Nf) x U(Nf) chiral flavor symmetry which is spontaneously broken. In order to 
realize the flavor symmetry, one introduces Nf D8 branes and Nf D8 branes transverse to 
X4 and localized [13]. In the probe approximation, where N c ^> Nf, the effect of the flavor 
branes on the geometry is negligible. The gravity approximation is valid when A5 ^> p. 
The dual theory in this limit is a 4+ 1 dimensional maximally supersymmetric theory with 
gauge group SU(N C ), compactified on a circle (with supersymmetry breaking, antiperiodic 
boundary conditions on the S 1 for the fermions) and coupled to Nf left handed fermions 
and Nf right handed fermions in the fundamental representation of SU(N C ). At energies 
lower than Mkk — 1/p the dual gauge theory is effectively four dimensional confining, 
with mass gap of order M^k- At the opposite limit A5 <C p, the theory is a 3 + 1 
dimensional and approaches a pure Yang-Mills theory where the mass gap is exponentially 
suppressed compared to 1/p. 

The D4 branes extend along x^,x^ while the D8 branes along the coordinates x M , f2 4 
and u. At the boundary the branes are localized at x 4 = (D8) and at rr 4 = L 8 (D8). 
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As they extend in the bulk they meet, and the chiral symmetry breaks since the group 
U(Nf) L x U(Nf) R which exists close to the boundary, breaks to the diagonal U(Nf) inside 
the bulk. By maximizing the flavor branes' distance to L$ = up at the boundary, the 
corresponding equations force them to meet at the lowest possible radial point u = Uk- 
Notice that in this setup the constituent mass of the quarks tends to zero. A gauge 
theory where the dynamics of glueballs and mesons influence each other can be obtained 
holographically by taking into account the backreaction of the flavor branes on the D4 
geometry. As already mentioned, there are several methods to consider this backreac- 
tion and take into accounts its effects. For our purposes we will focus on the localized 
backreacted flavor branes approach. 

To find the backreaction in the case under study one has to consider the D4-D8 
system in the Romans massive type IIA supergravity 4 . The corresponding equations for 
the dilaton and the graviton have been found [14], and as expected they contain delta 
functions. To solve these equations analytically in the full range of u is a very challenging 
problem. Instead of that, one can take advantage of the fact that the delta functions are 
codimension 1 and try to solve the equations perturbatively. The perturbative approach 
in this case seems to be consistent since the Green's functions are linear with respect to 
the distance of the brane and hence finite. This perturbative analysis close to Uk has been 
done in [14]. The small flavor expansion parameter qf used is defined as 

q f = g s K l0 K s = j-j- (6) 

where K 8 is proportional to the 8-dim brane tension and the Kio is proportional to the 
10-dim Newton gravitational constant. Notice that qf has dimension of inverse length. 
The solution when u — > is parametrized as 

ds 2 = e 2A l{ u, X4 ) (_ dt 2 + dx 2j + e 2A 2 (u, X4 ) dx 2 + JMw) du * + ^(u,*^ ( ? ) 

<p(u,x 4 ) = ^(f>(u,x 4 ) + 2A 1 (u,x 4 ) + 2At(w,x 4 ) , F [A) = Q C V A , 
where Q c = 3R 3 /g s = 3iiNJ^. The functions in the exponentials are expanded as 

Ai(u, x 4 ) = A U:i (u) + q f A dji (u, x 4 ) , (p(u, x 4 ) = (p u (u) + q f (frd(u, x 4 ) , 

where the u subscripted functions are the usual undeformed solutions of the g/ = 
equations, i.e. their exponents are equal to the metric elements of Eq. (3) and the index 
i runs from one to four. The functions determining the backreaction around Uk turn out 

4 Recently it has been argued that in the UV might be impossible to complete the WSS model. In this 
limit the massive IIA string theory between the D8-branes seems that can not be strongly coupled while 
remaining weakly curved [21]. 
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to be 




and the corresponding metric elements 




Using this background we examine the flavor effects in the static QQ potential in the next 
section. 

3 Static potential in presence of flavors 

As already mentioned, the full backreacted solution is not known and only the backreac- 
tion close to the tip of the cigar has been calculated. However, based on the fact that the 
string we are interested in lies almost entirely close to u k when u ~ u K , with u being 
the turning point of the string world-sheet, we claim that our results produce in this limit 
the leading corrections to the static potential. 

The open string between the two heavy quarks, is localized in the X4 dimension at 
an equal distance from the D8 and D8 branes. The string is attached to the boundary 
and hangs inside the bulk down to a minimum value u . In our analysis this minimum is 
chosen to be close to u k so that most of the string lies close to the tip of the geometry. It 
is then natural to think that the leading flavor effects for the particular WL come from 
the region of the world-sheet near u k and we indeed find that this is true. 

The flavor backreaction in our approximation should be perturbative and close to u k . 
This is the case close to u k , but as we depart from u k to higher values of u the elements of 
the backreacted metric (7) diverge relatively fast from the undeformed elements. Although 
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for our solution the important part of the metric is close to the tip of the cigar geometry, 
we would like to make sure that the string solution we find is valid in the deformed space. 
As we will see next, a localized string at x^ = vrp/2, solves the undeformed as well as the 
deformed equations of motion. Hence to summarize the picture of our string configuration, 
we consider a localized string at rc 4 = vrp/2 in the static gauge where it has a turning 
point at uq ~ Uk- It solves the equations of motion in the perturbative background as 
well as the original background, and it lies almost entirely close to the minimal point Uq 
where we argue that the flavor effects there will be the dominant ones. 

Along these lines we can perform also numerical computations. In numerics we would 
like to integrate through the whole background. So we need to introduce a 'cut-off' in 
the deformed metric to eliminate the flavor effects as the distance from Uk increases and 
we go outside of the validity of our approximation. One type of such cut-off could be the 
Heaviside step function, which would eliminate the effects in a non-continuous way. This 
could create a minor cusp in the world-sheet solution. Another way for the numerical 
integration is to use an exponential type 'cut-off' 

^ii=^ne 2 ^ lCxp[ - K{M - b)] , (11) 

and so on for the other elements (g u u is the undeformed metric element). The constant 
k defines how fast the flavor effects reduce as the distance from Uk gets bigger. The 
parameter b measures from where these flavor effects start to reduce. It is obvious that 
b = uq ~ Uk in our case. To determine the value of k, we have to set the (small) maximum 
distance from Uk within the regime where the perturbative solution is still valid and outside 
of which the flavor effects should disappear. It turns out that for the most of the cases 
k = 1 works well. We continue our numerical analysis in section 4. 



3.1 Flavor corrections for long strings 

Here we calculate analytically the leading flavor corrections in the static potential for long 
strings. Following the standard method we consider the string ansatz 

t — r, X\ — (7, x 4 = 7rp/2, u = u(a) , (12) 
where the equations of motion come from the Nambu-Goto action 

S = ~2~i J dcrdr^-g 00 (g u + g M x' 4 2 + g uu u' 2 ) , (13) 

and the metric elements are the ones of Eq. (7). The equation of motion for rr 4 is satisfied 
for rr 4 = vrp/2, so can proceed with the usual Hamiltonian formalism for the equation of 
motion of u and using the general relations derived in Appendix A. From Eq. (62) we 
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solve for the constant c which is the Hamiltonian and a constant of motion, in terms of 
the turning point u . The result reads 



U T ( If V u o - u l ( Au o + 5m 



R 3/2 



27u 



7/2 



(14) 



The integral for the string length including the correction due to the backreaction gives 



L — 2 I du 



^( 1 V"Q-"fc( 4 "0+ 5 "fc) + ( 191 "fc- 20 " 3 )y / " 3 -"fc) 



7/2 



R 3 u 3 



[u 3 - u 3 ) [e 



772" 



2 qfy /u3-ul(sul+4uZ) 



(15) 



Ur 



u 



where expanding in qf and keeping only 0(qf) terms we obtain 

L = L + Lf , 
with 

f°° I Rhi 3 
L °- 2 L dU \l(u 3 -u 3 )(u 3 -ul) 

being the corresponding integral obtained in the quenched version of the model and 



(16) 



(17) 



J f 



3/2 f°° 1 

Ju 135u 7 j/ 2 (u 3 — Mq) 3 ^ 2 a/ {u 3 — u\ 
\]u 3 -u 3 k (l91ugu| + 60m 6 - 20w 3 w 3 J - 141m 3 m| ) 

y/v*-u 3 k (-50u 3 ul-A0u 3 u 3 o ) 



+ 



the correction from the backreaction. 

Treating separately Lq and Lf, we first perform the expansion uq around Uk and keep 
until quadratic terms. Then we integrate the result and by keeping the leading terms at 
the limit u — > u k we obtain 



R 3 ' 2 (33 - 8VSn - 48 log (u - u k ) + 24 log (3u£)) 



72 v /m7; 

where we have also added the constant term 

R 3 ' 2 



zTT 



(19) 



(20) 



y/3uk 

which comes from the boundary. The leading contribution from the turning point is a 
logarithmically divergent term. Turning to the flavor term Lf we again expand Uq around 
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u k and obtain the indefinite integral. By keeping the leading terms in the definite integral 
we get 

_ 4q f B?/ 2 u k 

f ~ ~ o /q / = ' W 

9Vdy/u - Uk 

which again gives a divergent leading term when u — > u k . 

To find the energy in terms of L we need to invert the total length L given from 
Eq. (16) and solve for u . It turns out that the result contains the Lambert W function 
or product logarithm 5 



u = u k + 

108LambertW 

where 



. if u k , 
3 7/4 



2 ' 



(22) 



_ 11 7T ^ 3Ly/uj 

B - 16 2v/3 2i?3/2 • ^ 
Although the Lambert W function can not be expressed in terms of elementary functions, 
it can be Taylor expanded in qj giving 

M = u k (l + V3e B ^j ± 2 • 3- 5/4 q f u 2 k e% . (24) 

Using the well known fact that the constant c of Eq. (14) is the static force, we substitute 
the above 'plus' sign solution into Eq. (14) and integrate with respect to L obtaining the 
energy as 

2na'E = C + (|) ^ L - V3u k e B - ^qfU^ , (25) 

with C a constant. From the above relation it is obvious that the correct solution in 
Eq. (24) is the minus sign solution, since this gives the physically expected predictions 
and it is in agreement with the numerical results. Moreover the plus sign expression of 
the Eq. (24) gives contradicting L(u ) behavior to the Eq. (16). 

The 'minus' sign solution does not give any correction in the leading order of qf. Hence 
the energy reads 

2W£ = C + (^) 3/2 L - y/3 Uk e B + O {qje cB ) , (26) 

where c is a constant and the qj correction term should come with a positive sign in 
order for the static force to be screened. We see that the flavor correction to the energy 
is similar in form to the correction present already in the undeformed background: they 
are both ~ e~ L , although they have different coefficients. Evidently, at this level of 
approximations, we do not see the expected from lattice QCD flavor correction ~ 1/L. 



5 The Lambert W function is a set of functions which defined by solving the equation z = W(z)e w ^ z \ 
for any function W. This function is also called product logarithm and can not be expressed in terms of 
elementary functions. 
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3.2 Flavor corrections for shorter strings 



3.2.1 Probe approximation limit 



We would like to examine also strings of physical length of order one. This could be 
relevant in the case where the string breaks due to a vanishing pion mass at a regime 
where the static energy is around zero. In this section for convenience we define y = u/uq 
and A = Uk/uQ. In the absence of dynamical flavors we have to solve the following two 
coupled equations: 

£==^ = ^f , dy (2-) 



and 



u k A 



dy 



- fdy- 1 



(2? 



In this dimensionless basis, the infinitely long string (A — >■ 1), has tension = 2. 
For A < 1 we can expand as 



with 



{ A^ 

y 3 



Ck 



( y4. 



3\ k 



-1/2 
k 



(29) 



(30) 



the binomial coefficients. Then, using the equation 



dy 



,3+A 



vW - 1) 



— my 



E 

fc=0 



3\ k 



k=0 



m 



y/nF(k- 1/3- A/3) 
A + 1-3A; ' 3 I\fc + l/6-A/3) 



+ 



we can express the two integrals of L and E by setting A = —3, m = and A = 0, m = 1 
respectively, as 

f _ v 7 ^ £(fc + 2/3) 3fc+1/2 

- 3 2-^ + 7/6)^ ' 
~ 2v^F^ r(ib-l/3). 



3fc-l 



(31) 



From now on we drop the ^-independent constant from E and by defining 



Ofc :— Cfe- 



0r T(k + 2/3) 
3 r(Jfe + 7/6) ' 



h ■= c k 



2y/nT(k-l/3) 
3 r(fc + l/6)' 



x 



:=VA, (32) 
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we can rewrite the system as 



k=0 

oo 



E = J2 b ^ k ~ 2 - ( 33 ) 



fc=0 



Taking derivatives of both Eq. (33) with respect to x it is easy to show that 

n ,dE /M \3/2 

a relationship which actually holds more generally (see for example [23]). 
Next we have to invert the series of L. We first compute the derivative 



dx 



(35) 



dL Em^ + I)^ 
The solution for x can be then written in terms of an expansion around x = x c 

x = d + dj + d 2 l 2 + dj 3 H , I := L - L c , (36) 

where L c is the value of L at x = x c , that is L c = YlT=o x t h+1(1 k- This immediately implies 
that do = x c . The n'th coefficient d n is computed by taking d n x/dl n and setting x = x c : 

i , . = i &M±M^ 

1 ' ZT =0 (M + i)xfa k ' d2 ■ 2 CLT=o( 6k + 1)*^) 3 

and so forth. Substituting this solution in the series for E one obtains 

E = - b ° +bl ( Xc + dl i+d 2 P + . ■■)* + ■■■ (38) 

(x c + dj + d 2 l 2 + • • • ) 

and putting back the units one obtains the expression for E(L). An observation is that if 
we expand the above expression in small /, the coefficient of the term linear in / converges 
to the infinite string value a*- -' = 2 from below: by summing 10 5 terms in the series, the 
coefficient is 1.999772328. 

An interesting first application of the above formulae is when x c = 0. This corresponds 
to strings of nearly zero length, i.e. to the UV of the gauge theory, that is to its 5-d 
Coulomb phase. Then it is easy to see that the leading contribution to the energy is 
E = which means that 

n : " 2 / r(2/3) \ 3 1 
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or in terms of the dimensionful quantities that 

4 /^r(2/3)\ 3 i? 3 



E 



27 V r(7/6) J 2ira' L 2 
2 f ^T(2/3) \ 3 glN c 1 

27 V r(7/6) ; 4tt 2 L 2 ' 1 J 

The coefficient d 5 of the five-dimensional Coulomb potential 

<«) 

was computed in the mean-field expansion of a five-dimensional pure SU(2) lattice gauge 
theory, near its bulk phase transition [24] . The lattice coupling, which is the only dimen- 
sionless parameter on an infinite lattice, is defined as (3 = where a is the lattice 

spacing. The difference between this theory and the unflavored Witten model in the UV, 
is the adjoint scalar field that the latter has, but this is not expected to affect very strongly 
the Wilson loop, which is an object constructed out of gauge fields. The result of the 
lattice computation was d 5 /a = 0.0626 for /3 ~ 2 [24]. Let us then consider the Witten 
static potential in the UV Eq. (40) (see also [22]) and identify its Yang-Mills coupling g$ 
with the corresponding lattice coupling g 2 N c = ^pa. This should be a valid approxima- 
tion as effects of finite lattice spacing become small and could allow us to make a first, 
crude but interesting comparison between two very different computations of analogous 
quantities. This leads us to the relation 

ya = ±(WW)'l-£ =0M6 «L. (42) 
1 27 V r(7/6) J 4tt 2 p p K ' 

For TV 2 ~ P there is an impressive agreement with the prediction of the lattice model. 
For large N c however this requires a large P that corresponds to the perturbative regime 
of the five-dimensional gauge theory, not near the bulk phase transition where P is rather 
small (and where the 577(2) lattice result comes from). Clearly at the moment no definite 
quantitative conclusions can be made but the comparison certainly motivates a gener- 
alization of the mean-field lattice computation to SU(N C ) for general N c , where a more 
direct comparison can be made. 

In order to probe strings of finite length however, we need to understand x c — 1, its 
maximum possible value. This corresponds to expanding around u = -u fc . Notice that in 
this case 

^ = —1:477^^- 0.70125 (43) 
3 f^ Q T{k + 7/6) v ; 

approximately. This corresponds to a string of physical length L c ~ 1.14. This is approx- 
imately also the length where the energy is zero, so with massless pions in the theory it 
should be around the string breaking scale. This could be an interesting point to expand 
around when we add flavors, which we do next. 
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3.2.2 Unquenched flavor 



By introducing flavors we write Eq. (21) as 



9^ 



c k x 



2k+l 



with e f := q f u k , 



(44) 



fc=0 



where ej is the dimensionless flavor expansion parameter. The previous probe relation- 
ships generalize to 



fc=0 



d 2 



£L>(6* + l)x?a fc - ^ Er=o(2^ + l)**c fc ' 

9^ 



(45) 
(46) 



Er=o(6A: + l)x?a* - ^ EZo^ + l)xfc k ) 



The maximal L c does not change significantly when e/ is small. The energy can be 
computed from 



1 ul - u\ , 2 («g - «D 3/2 



27 



u 



7/2 



(47) 



We first rewrite it in terms of dimensionless quantities as 



dE _ 2 
dL ~ x 3 



1 ( 1 



1 + 2£ 'l6 U^" 1 



1/2 



+ 



1 

27 



- 1 



3/2 N 



(48) 



and then we integrate term by term. The first term gives 

E = . (49) 

d l {x c + dil + ---) 2 

We kept x c general so that we can check against Eq. (39). Indeed it is easy to see that 
for x c = the results agree. We define d w and d±f 

1 



dw + e/di/ 



(50) 



which can be identified in a straightforward way from Eq. (46). Then, performing also the 
integrals proportional to e/ in Eq. (48) and expanding in e/, to leading order we obtain 



E = 



d 3 
u io 



{x c d 10 + I) 2 



d\ Q du 



+ 



2,d 2 Qdnl 



(x c d 10 + l) 2 (x c d 10 + l) 3 
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+ 



(51) 



where we can now set x c ~ 1. If we expand in small / and sum a large number of terms, 
we find a string tension whose unflavored value is corrected only by small flavor effects. 

Finally we would like to have an improved analytical handle on strings, with 0(1) < 
L < oo so we include a higher order correction to the long string expansion of Eq. (16): 

L (2) = _±( 13 (-53 + 42V3n) ^_ 13 log(s) \ 
f 2 y 204120 108^ 1620 J 1 1 

where 

a parameter that makes sense for < z < oo. The length around which we are expanding 
the string is now 

L c = a 2 log z c - a 3 z c + (3iy/zZ (54) 

with 

1 13e 2 , e) , 2e f 

d2 = 3" 3240' d3= 2W!' * = -9^ (55) 
Clearly, since only in the leading terms we have kept uq ^ Uk, one should think of z c being 
large. By performing similar steps as above, when / is small we obtain 

3^ 9 ,2-S 3 / 4 20 -3 1 / 4 , 

E ^ 2+ — + i^<vr + ^ )+ '"> ;+ "' (56) 

where I — L — L c . Using this expansion L c can be lowered significantly but cannot be 
taken too low because in order to derive this result, we have set u = in sub leading 
terms. We can see that as z c increases as we approach the infinite string tension value 2 
(with small flavor corrections) from above. 



4 Numerical analysis 

The numerical analysis can be done either by using the method with the cut-off of Eq. (11) 
discussed in the previous section or by introducing a step function which defines the region 
within the backreacted flavor metric is valid. However, one should be careful with the 
values of the step function and especially of the cut-off otherwise the integrals might not 
converge. Both of the methods give similar results. 

By doing the numerical integration for both E and L and using the renormalization 
scheme of the subtraction of the infinite mass quark we obtain the results which are 
presented in Fig. 1. 
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Figure 1: The energy of the QQ system with numerical integrations in the flavored backreacted 
background and in the original one. The flavor effects in our analysis come as perturbations so 
the results are expected to be close at least for large L which correspond to world-sheet turning 
points close to the u^. We observe that close to Uk the undeformed energy is slightly less than 
the deformed one. The options we follow for all the plots are: The values that are chosen 
are qjUk = 0.15, uq/u^ ~ 1 + 5 • 10~ 5 and increasing slowly close to this region. The flavor 
backreacted metric is valid only around Uk- The energy and the length are divided everywhere 
by Uk- With the red dashed line are the results of the quenched WSS model and with the blue 
solid line the results of the backreacted WSS model. 

The static force between the QQ is plotted in Fig. 2. 6 On the left we show the result 
for the WSS model for Nf = and Nf ^ 0. On the right we show the result of the lattice 
QCD analysis for Nf = and Nf = 2. The Nf = lattice data is taken from [26] and 
the Nf = 2 data from [3] and [27]. In the two sets of data the extrapolated to zero quark 
mass value Lo/a = 7.05 has been kept fixed, with Lq the Sommer scale and a the lattice 
spacing. In the WSS model we have also kept a scale L fixed, however it is not possible to 
unambiguously determine its value. In the plot, we have chosen the crossing point of the 
flavored and non-flavored data of the WSS model to coincide with the one of the lattice 
data. There is an obvious qualitative similarity between the two plots but at this stage 
it is not clear if the fixed scales in the two approaches have the same physical meaning. 
For large L/Lq, where we know the backreaction effects, the static force in the flavored 
background is less than the one in the non-flavored case. This could be due to the fact 
that in the presence of flavors, there are screening effects caused by the virtual quark-pairs 
produced between the heavy quarks, and reduce the color force. Our interpretation of the 
plot for smaller values of L/L relies on the validity of the backreacted solution in that 
regime. 

It is interesting to observe that the strings with the same length L, have lower energy in 
the undeformed background than in the deformed one at least in the region where w ~ u k- 
Moreover, string world-sheets that are equally radially extended from the boundary have 

6 Alternatively we could have plotted the renormalized charge agg(L) = L 2 F(L)/(4/3). 
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Figure 2: The static QQ force, as a function of the heavy quark distance. On the left the 
results obtained from our holographic analysis. On the right the lattice results for Nf = and 
Nf = 2 taken from [26], [3] and [27]. With the red dashed line (red points) are the results of 
the non-flavored case and with the blue solid line (blue points) the results of the flavored case 
of the WSS model (lattice analysis). In the lattice data Lq is the Sommer scale. In the WSS 
data it is some fixed scale. 

lower energy in the flavored background, and smaller lengths L are required to reach the 
minimum distance uq. Hence, to obtain the same QQ energy in the flavored and non- 
flavored background, the string in the flavored case need to be extended less and have 
smaller length than the one in the original background (Fig. 3). We expect that these 
observations carry on in the case of the full backreacted solution of the background. 

We observe that since our flavor effects are perturbative the modifications to the 
backreacted Wilson loop are small. However there are important differences as the turning 
point goes far away from Uk and this happens because there our approximation stops to 
be valid and because of the mass subtraction term which still gets contributions from the 
region close to Uk- When uq gets outside the region where the flavor effects are taken into 
account, the first terms of the equation Eq. (65) do not receive any flavor contribution, but 
the last one coming from the mass of the infinite quark does get important contributions no 
matter how far from Uk the uq is. This situation would be different, if for the regularization 
scheme the Legendre transform was used. However to use this renormalization to subtract 
the divergence, the conditions for the background derived in [28] must be satisfied. In 
our approximation it seems that these are satisfied at least approximately but would be 
interesting if one could check it for the whole backreacted background. 

5 Discussions and Conclusion 

We have touched upon several problems that obstruct the WSS model from describing 
QCD. We have focused on the fundamental properties encrypted in the static potential and 
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Figure 3: On the left the energy of the string with respect to minimum point uq. On the right 
the length dependence of the string with respect to no- The string worldsheet in the flavored 
background reaches the same energy with the string in the undeformed one for lower values of uq. 
Furthermore, the string in the flavored background that extended to a minimum point uq has 
length L smaller than the string with the same minimum point in the undeformed background. 
With the red dashed line are the results of the WSS model in the probe approximation and with 
the blue solid line the results of the backreacted WSS model. 

static force rather on shortcomings at the phenomenological level. We have seen that the 
regime where the flavor correction to the Liischer coefficient develops that the lattice seems 
to observe between the Sommer scale and the string breaking scale, cannot be located. 
We consider this particular flavor effect in lattice QCD special because if it admits a string 
description, it is an order 2 effect in the Liischer coefficient. For that reason, we have 
scanned analytically via two different expansions and also numerically the WSS Wilson 
Loop at essentially all scales. We have seen no sign of a regime where flavor corrections to 
the string tension in inverse powers of L develop. We have argued that a' corrections are 
unlikely to change this picture. According to this discussion the perturbative effects on 
the backreacted background may not be enough, and possibly a full backreacted solution 
is needed in order to understand via holography the phenomena observed in lattice QCD. 
In case these problems persist even if a fully backreacted background is found, it would 
mean that a more drastic modification of the WSS model is necessary. 

The color force beyond a certain scale reduces in the presence of flavors due to screening 
effects caused by the virtual quark-pairs produced between the heavy quarks. We have 
found this behavior for large distances between the quarks, where the validity of our 
approximation is best. This has been also noticed in the backreacted D3/D7 model where 
the screening effects result to a change of slope in the Regge trajectories [29]. A naive 
comparison between lattice data and our results shows qualitative agreement as depicted 
in Fig.2. Nevertheless, at this point one should be aware of two sources of uncertainty. 
One is the validity of the range of the backreacted solution which could result in the 
low L/L part of the static potential not to be trusted. The other is related to the 
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identification of L as the Sommer scale in the WSS model. 

A full understanding of screening in the WSS model and more generally of its physical 
scales is subtle. The origin of the subtlety is an obstruction in the precise determination 
of a physical length scale: in the unflavored model the absence of a Sommer scale and 
in the flavored case in addition a difficulty in determining the string breaking scale. The 
former obstruction is due to the fact that as the theory is followed from the IR towards 
the UV one does not seem to pass through a regime where it is 4-d and asymptotically 
free, nor crosses a first order phase transition. Instead it passes smoothly into a 5-d 
Coulomb phase, as Eq. (40) indicates. There is not much to say about this problem as 
at present seems to be an intrinsic property of the WSS model. We could only suggest a 
possible way out which involves the construction of a background dual to an anisotropic 
5-d gauge theory as opposed to an isotropic one, as it is the case here. The discussion of 
the latter obstruction is also subtle. To begin, for the same reason as in the unflavored 
case there is no Sommer scale. One would hope that the backreacted background contains 
information about the breaking of the string, which would give us an alternative way to 
extract a physical scale. In the Veneziano limit of the MN background 7 system where 
backreaction of the smeared flavor branes have been considered [30], it has been noticed 
that the backreacted background contains the information of the string breaking and as 
the number of flavor branes is reduced compared to color branes, the length where the 
string breaking happens increases for massive dynamical quarks. Hence, one could think 
that in our setup, because the backreaction comes only as perturbation to the original 
background, string breaking happens at very large L close to infinity, where due to the 
constraints in numerical integrations we could not see it in the plots. In our setup however 
the D8 and D8 branes are placed antidiametrically on the circle. Therefore the constituent 
quark mass is zero since it is given by the energy of a string extended from the minimal 
position of the background to the minimal radial position of the corresponding flavor 
branes. Moreover, the pions in our model are massless and one would rather physically 
expect that string breaking happens at short lengths, since a very small amount of energy 
is needed to create pions. According to this argument and our results, the string would 
break in this case around L ~ 1. Hence there seem to be two opposite ranges on the scale 
of string breaking and in our approximation it seems not to be possible to get a conclusive 
answer just from our model. 

As a result, the length scale around which one should expand could possibly be far 
from the domain of validity of the backreacted background. This is not necessarily the 
case of course, but it is not possible to perform a reliable quantitative treatment until a 
better understood version of the model with massive pions is constructed. In principle it 

7 The analysis in the D3/D7 system can be found at [31]. 
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is possible to incorporate some mechanism that gives vev for the chiral condensate and 
the quark mass [32]. There are two main difficulties that we see at present with this. One 
is that we do not know of a formalism fully under control for Nf > 1 and the second is 
that one should take into account the resulting modifications in the backreaction. 

Concerning the behavior of the string world-sheets we have found that the strings with 
the same length L, have lower energy in the undeformed background than in the deformed 
one at least in the region where uq ~ u\~. Moreover, string world-sheets that are equally 
radially extended from the boundary have lower energy in the flavored background, and 
need for smaller lengths L to reach the minimum distance uq. As a result to have same 
QQ energy in the flavored and non flavored background, the flavored string is extended 
less and have smaller length than the one in the original background. These observations 
should carry on in the full backreacted solution too. 

A possible interesting analysis of the static potential could be done also away of Uk, 
for short string which stay close to the boundary. There the perturbative backreaction is 
known [14], but the analysis seems to be doable only at numerical level. The backreaction 
at this range should be also similar to the NJL dual gravity theory [33], and that is 
an additional interesting point. However, the most challenging problem would be to 
obtain the full backreacted solution of the WSS model and even better by including the 
bifundamental tachyon field in order to have massive pions. This could be done either by 
considering the localized branes or by the smearing techniques. If the full solutions were 
obtained, for small backreaction (small number of flavored branes compared to the color 
ones) our results should be reproduced. Additionally, in this background one should be 
able to study explicitly the string breaking and make a comparison with the lattice data. 
The extension to the finite temperature backreacted background is also of big importance 
especially for quark-gluon plasma physics since the effect of the dynamical quarks to the 
plasma physics should be taken into account. 
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Appendix A: Strings in general weakly coupled backgrounds 

In this section we derive the turning point of world-sheet of a string in static gauge 
in weakly coupled backgrounds and their energy. The confining string in the QCD back- 
grounds, is a special case of the strings we examine here. The following equations have 
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been derived in various forms in several in other papers too [15, 18]. By writing the metric 
of the space as 

ds 2 = goodr 2 + gudx 2 + g uu du 2 (57) 
we choose the static gauge for the string 

xq = t and x\ = a, (58) 



which is extended in the radial direction too, so u = u(a). Then the induced metric 

Goo = 9oo, G u = g u + g uu u 2 . (59) 



G a /3 = gMNd a X M dpX N for our world-sheet read: 



Supposing that we are working in Lorentzian signature the Nambu-Goto action is 8 

S = 7^7 J dadT V / -g 00 (g 11 + g uu u' 2 ) =: J dadry/D . (60) 
The Hamiltonian then is equal to 



H= googu (61) 



and is a constant of motion. Setting it equal to c we can solve for u' and get the turning 
point equation 

V c 2 g uu 

which is solved for 

9uu = , or 9n = , or g 00 g u = -c 2 . (63) 

The above equations, normally the last one, specify how deep the world-sheet goes into 
the bulk, and we call the value of the turning uq. 

The length 9 of the two endpoints of the string on the brane is given by 



, , f U ° du _ 2 r du / ~ 9uu ° 2 (64) 
u' Ju V (doogu + c 2 )g u ' 



8 In the case of Euclidean signature, the formulas change with a minus sign wherever the 500 element 

is. 

9 The limits of the length integral depend on where we choose as the starting point measuring L. 
When the string in the boundary extends from — L/2 to L/2, then the corresponding solution of the v! is 
positive for (0,L/2) since the turning point corresponds to L = 0. When the string extends from (0,L) 
the u' in (0, L/2) is negative. In any case the final result in the definite integral is Eq. (64). 
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Moreover, the energy of the string using as renormalization method the mass subtraction 
of the two free quarks is 



2na'E 



f oo /*oo 
'uo J UQ 



= cL + 2 



du^-g uu g 00 



1 + 



uo 



9ii9oo 



du^-googu 



.(65) 



Notice that we already used that fact that the world-sheet is symmetric with respect to 
turning point U and hence the RHS of the above equation are already multiplied by two. 

Appendix B: The D4 Wilson loop in the IR 

In the static gauge X = a, X 1 = r the metric Eq. (3) can be brought in the form 



ds 2 



1 



2T 



Y ■ Y) )dX ■ dX + dY ■ dY + R 2 dtti 



with dX-dX = dX\ + dXl Y-Y = F 4 2 + F 5 2 , T = 2na'a^ and where 

1 

m b = - 



(66) 



(67) 



is essentially the "Kaluza-Klein" scale governing the lightest glueball. This immediately 
suggests us that only the transverse X 2 ^ fluctuations are massless and thus that the 
induced Liischer term in the static potential is likely to be the same as that of a four- 
dimensional non-critical bosonic string in a flat background, cf augc = —n/12. Next, we 
would like to actually perform the integral 



e -s eS [x] = / dYe - s ^ . 



(68) 



The result is [19] 



l + -(d a X).(d a X) + 



S cS [X] = <r (1) / d 2 a 
Jm 

+ Va^) [ da°[pL B + b 1 (d 1 X)-(d 1 X) + ---} 

JdM 



(69) 



where the dots stand for the world-sheet fermions and higher derivative terms and 

crW — cr^(l + AT) (70) 
is the corrected string tension renormalized by 



AT = — 

8tt 



/ 



(71) 
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The m/ are world-sheet fermion masses. There is a wave function renormalization involved 
to bring the result in the above form, which is finite provided that the condition X^6 m b = 
^2 j rrif holds. This condition holds for most known holographic backgrounds, including 
the one under study. The constant can also be computed 



1 



77V 



(72) 



where the index a counts Dirichlet and a' Neumann boundary conditions. Finally, the 
coefficient in the boundary terms vanishes: b\ = 0. All the above implies that after 
integrating out also the X-fluctuations, the static potential of the Witten (and in fact of 
a fairly large number of holographic) model in the IR can be summarized as 

E m = const. + a^L - ^ + 0(L~ 3 ) . (73) 

Evidently while the string tension is renormalized, the Luscher coefficient retains its uni- 
versality. The coefficient — tt/12 assumes that the perturbatively zero modes on the S 4 [34] 
pick up a mass at the non-perturbative level [19]. 

The perturbative backreaction used in this paper, if it has any effect on the Lscher 
term it will probably give a small mass to the massless transverse bosonic fluctuations and 
perhaps correct by a small amount the already massive fields. Thus, the corresponding 
term could change into a term of the Yukawa form 

(-7T/12)e- mi 

L 

with 777 the small mass due to the backreaction. Clearly, this will be a small effect. 10 



(74) 
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